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. We investigate the entanglement content of the ground state of a system characterized by effective 

T-H ' elementary degrees of freedom with fractional statistics. To this end, we explicitly construct the 

\ ground state for a chain of A'' spins with inverse square interaction (the Haldane-Shastry model) 

. in the presence of an external uniform magnetic field. For such a system at zero temperature, we 

5_( ' evaluate the entanglement in the ground state both at finite size and in the thermodynamic limit. 

^ I We relate the behavior of the quantum correlations with the spinon condensation phenomenon 

, occurring at the saturation field. 

PsJ ■ PACS numbers: 11.25.Hf, 74. 81. Fa, 03.65.Yz 

CD ; I. INTRODUCTION 

i: 

C/5 , Quantum states of interacting many body systems are inherently entangled. This simple statement has motivated 
' a wide physics community to study, in the light of quantum information theory, models typically explored in quantum 
\ statistical mechanics Besides the motivation of the crucial role played by the entanglement as a resource for 
quantum computation protocols, the promise is to acquire a deeper understanding of the physical systems themselves, 
with possible spin off on open problems in condensed matter. There are indications that, while local interactions are 
possibly the most common mechanisms leading to non local correlations, also statistics can produce entanglement. 
Q ' Indeed; because of the symmetrization constraints, most gencrically identical particles are not in a product state. 
^ , Even though the resulting many body wave function cannot lead to genuine entangled state because of the lack of the 
tensor product structure in Fock spaces 0], Vedral Q demonstrated that a free fermion gas enjoys spin entanglement 
on distances of the order of the inverse Fermi momentum; on the contrary, no such kind of entanglement has been 
. found for the polarizations of identical bosons. An alternative route to capture non local correlations in the states of 
identical particles was pursued by Zanardi and coworkers through the so called mode- entanglement, directly defined 
in the Fock space Q. As a basic support to this point of view, one can argue that separable states in the dual space 
■ correspond to configurationally entangled states. Furthermore, in favor of such 'statistics entanglement', it should be 
, noticed that several protocols have been suggested to exploit it for computational purposes [H . 

' In this paper, we deal with statistics in a one dimensional interacting many particle system. In one dimension, 

' statistics is peculiar as particles need to scatter each other in order to exchange their position. Accordingly, one 
dimensional statistics can be formulated as a boundary condition for the many-body wave function in the configu- 

' rational space @. In this context, particularly relevant are the Calogcro-Suthcrland models (CSM)s, i.e., quantum 
^ . mechanical models in one dimension where the interaction among particles is inversely proportional to the square 

' of their distance. Indeed, such models can be regarded as describing a gas of identical free particles with fractional 
1^ , statistics 0, [3i which has made them relevant in the understanding of the fractional quantum Hall effect. 

' The main motivation of our work is to investigate the kind of entanglement emerging from fractional statistics. 
We note that, because of the lack of any tensor product structure, the entanglement between particles with fractional 
statistics has not been fully quantified yet[^. We concentrate on the spin entanglement emerging from quasi-particles 
with fractional statistics. To rely on some tensor product structures of the underlying Hilbert space hosting the 
quantum system, we consider the Haldane-Shastry [loj model (HSM): an anti-ferromagnetic spin-1/2 chain with 
inverse square exchange interaction among the spins. 

Also the HSM (in the dual space) can be interpreted as an ideal gas with fractional statistics [lH . The connection 
with the CSM can be put on firm grounds following the scheme developed by Polychronakos; namely, first introducing 
some internal degrees of freedom into the CS model and then employing the so called 'freezing trick', to fix the 
particle positions [l^l- As can be understood by a direct inspection of the eigenstates of the two models, it is 
worth mentioning that the correspondence between the CSM and the HSM holds only for a well-defined value of the 
parameter characterizing the statistics of the excitations in the CSM, identifying a 'scmion' gas (i.e., a gas of particles 
with statistics half than the Fermi statistics) . One of the most spectacular evidences that the HSM is indeed an ideal 
gas of semions is provided by the response of the system to a local perturbation: it can be proved that a single spin flip 
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causes the creation of a pair of collective excitations (spinons). Higher 'angular momenta' excitations are generated 
by two or more spin flips [l3| . 

In the present paper, we study the entanglement in the ground state of the HSM in an external magnetic field. 
In fact, the applied magnetic field causes a certain number of spin flips to occur in the ground state: this allows us 
to study quantum correlations in the presence of spinon condensation. We shall demonstrate that the divergence of 
the spinon 'Fermi' wavelength occurring when all the spinons in the system actually condense (saturation field), a 
characteristic trait of the semionic statistics, causes a divergence of the entanglement length. We also note that in our 
work wc study the interplay between the range of local interaction and the range of the entanglement. We eventually 
find that a specific pattern of entanglement characterizes the system with a finite range interaction. 

The paper is organized as follows: in the following section, we introduce the model and explicitly construct its 
ground state as function of the magnetic field; in Sect. IIIIl we employ the single- and two-spin ground state averages, 
computed in Appendix |X] to evaluate various entanglement measures. Finally, Sect. IIVI contains some concluding 
remarks. 



II. THE MODEL 

The Haldane-Shastry (HS) Hamiltonian describes a one-dimensional spin-I/2 chain wrapped around a circle of 
unit radius, with an anti-ferromagnetic interaction whose strength is inversely proportional to the square of the chord 
between the corresponding sites. When the model is placed in an external uniform magnetic field h, its Hamiltonian 
is given by 

In Eq. ([T]), Sq, is the spin operator and the spin coordinate of the site a = 0, l,...,iV— 1, with N the total number 
of spins. Periodically boundary conditions = Za^^, and S„ = Sq+at, allow to parameterize the positions on the 
circle by the N^^ roots of the unity Za = e^'^*"/^. For h = 0, the ground state of the system is a non-degenerate spin 
liquid with total spin (a spin sing let) [3. At zero field, flipping a spin in the ground state corresponds to creating 
a (fully polarized) pair of spin-I/2 collective excitations [3], dubbed "spinons" by Haldane 0. The spinons keep 
their integrity when scattered off each other, so that they can be thought of as true quantum particles, despite their 
collective nature. Any time one more spin is fiippcd, an additional spinon pair is created in the state of the system; 
hipping K/2 spins is, thus, equivalent to creating a fully polarized A'-spinon excited state. The Zeeman energy in 
Eq. ([T]) lowers the energy of the fully polarized states, since h works like a chemical potential for the spinons. Thus, 
upon increasing the magnetic field the (even) number of spinons K increases and one of the fully polarized A'-spinon 
states becomes the ground state of the system. This state can be written as 

|$gnd,A')= <^>k{zi,...,zm)\zi,...,zm) , (2) 

where zi, . . . , zm denote the positions of the f-spins in the state, the remaining ones pointing downwards. In Eq. ([2|), 
the number of f-spins is given by M = {N — K) /2, where AT is the number of condensed 4,-spinons in the ground state; 
the state |zi, . . . , zm) is understood to be equal to |zi, . . . , zm) = . . . | \.)^^ and the amplitude ^xizi, . . . , zm] 
has the form 



^k{zi,...,zm) 



The number of condensed spinons is determined by the value of the external field. In order to find out how AT 
increases with h, one has first to calculate the energy Egnd[N,K,h] of the state |$gnd,-f'^), as a function of K and 
h at fixed N, and then to determine the value of AT that minimizes Egnd[N, K, h], at fixed h. Egnd[N, K,h] can be 
obtained, for example, by a straightforward generalization of the approach developed in Ref. [isj : since (zi, . . . , zm) 
is a homogeneous polynomial of zi, . . . , zm-, the summations over spin indices, contained in the expectation value 
{^gnd, K\H\^gnd, K) , are replaced by derivative operators that are understood to act onto the analytic extension of 
^fcizi-, . . . , Zm), hi which the spin coordinates are allowed to take any value on the unit circle. After computing the 
derivatives, the variables zi, . . . , zm are constrained back to the lattice sites. The ground state energy is then the 
minimum energy of a AT-spinon condensate at finite h; it can be expressed as 

Eg^4N,K,h] ^ E[qi{M), . . . ,qKiM),qK^,iO), . . .QKiO)] ^ , (4) 
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is the energy part associated to the bare HS Hamihonian, and 
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(5) 



(6) 



denote the spinon (pseudo) momenta with < mi < . . . < mx < M. 

According to Eqs. ([I])-®, the minimmTi energy if -spinon fully polarized state is realized for mi = . . . = mj^/2 = M 
and mj^^2+i = ■ • • = mx = 0, which corresponds to the condensing spinons being equally distributed at the ends of 
the single-spinon Brillouin zone, as depicted in FiglT] When all the spins in \^gad,K) are polarized, that is, when 
K = N, the ground state of the system is given by the simple factorized state |$gnd,iV) = | i)®^. In the language 
of spinons, |$gnd,-^) corresponds to a state with the single-spinon Brillouin zone completely filled. In fact, the 
"saturated state" |$gnd,-^) becomes the ground state of the system as soon as h > hg, where hg is the "saturation 
field" (to be estimated below). The interval < h < hg is, therefore, divided into N/2 regions labelled by the number 
of condensing spinons. The n*^ region, with 2n spinons (n = 0, . . . , N/2 — 1), ends with the energy crossing between 
Egad[N, K = 2n, h] and E'gndl-^, K = 2n + 2,h]. The width of each interval becomes smaller, and the set of crossing 
points becomes denser, as the number of spins increases, which gives rise to a quantum instability similar to those 
found for the Dicke and the XX models [1^1 . 

Let K[N, h] be the number of condensing spinons, ob- 
tained with the minimization procedure outlined above. 
The variation of K[N,h] with respect to h/J for dif- 
ferent values of TV is reported in Fig. [2]A.. where it is 
shown that, for large enough A^. the steps in K[N,h] 
are smoothed down. This implies that the spinon den- 
sity K[N, h]/N, shown in Fig.[2P, can be approximated 
with a continuum function pk{N, h), whose dependence 
on h at fixed iV is given by 



Pk{N, h)^l 
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It is worthwhile noticing that px [N, h) is independent 
of A^, as iV — >■ cx), which is also displayed in Fig. [5)B. 
The value of the saturation field hg can be derived from 
Eq. ([7|), by requiring that px{N,h) = 1, for h > hg. 
As a result, one obtains hg sa tt^J (-j + 3^), which, for 
large enough A^ (> 10), gives hg w tt^ J/4. Accordingly, 
the number of condensed spinons as a function of the 
external field is given by 




FIG. 1: Single-spinon level occupancy corresponding to the 
state |$gnd, -K^). The spinons symmetrically populate low-energy 
states around the corners of the Brillouin zone (which ranges 
from — 7r/2 to 7r/2). The corresponding "spinon Fermi surface" 
is composed by two disconnected 'Fermi' points (the largest val- 
ues of the momenta at which we can create a spin pair), at 
Mo = {!-§). 



K[N, h] 
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for h < h^ 



(8) 



A^ for /i > hg 

This can be regarded as a "thermodynamic limit formula", which is confirmed by the plots of Fig. [2]B. 



III. GROUND STATE ENTANGLEMENT 



This section is concerned with the evaluation of two kinds of entanglement measures, both depending on the external 
magnetic field and on the number of spins. One of these is the one-tangle Ti{N,h), which quantifies the amount of 
entanglement shared by each spin with the rest of the system; the other is the concurrence of any pair of spins 
C{N, h, r) that also depends on the relative distance between the spin sites r = a — a' {a ^ a' = 0, . . . , A^ — 1), 
because of the translational invariance of the system (see also [l3|)- 
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FIG. 2: A: Number of condensed spinons, K[N,h], vs. h/J, for different values of A'^ ranging from 10 to 400. B: Density of 
condensed spinons p[N, h], vs. h/J, for A*' = 20 and — 100; the thermodynamic expression ((Tjl, holding exactly for A'^ — >■ oo, 
is also shown for comparison. 

Using the concurrence, the two-tangle T2{N,h) — ^^.C^iN ,h,r) and then the entanglement ratio R{N,h) ~ 
T2 (N, h) /ti [N, h) will be evaluated. They will provide information about the fraction of entanglement shared by 
pairs with respect to the total (i.e. both bipartite and multipartite) entanglement to which each spin participates. In 
order to evaluate these entanglement measures, the mean values of the spin operators and of the two-spin correlation 
functions are needed on the state j^'gnd,^)- 

The technical steps required to obtain these quantities are reported in the appendix where their explicit ex- 
pressions are given, and where the thermodynamic limit is thoroughly discussed. In particular, all of the one- and 
two-spin properties can be determined analytically for any value of the external field h, and for any given number 
of spins iV, except for the longitudinal correlation function. The latter is obtained numerically at specific values of 
h (see appendix 1X1 for the details). 



A. One-tangle 



To compute the one-tangle, one needs to obtain the average 
value of the three cartesian cornponents of each spin 5^, with 
j = x,y,z, on the ground state Since the applied magnetic 
field does not break the rotational symmetry around the z- 
axis, one gets 



(*g„d,if|5^|$g„d,^^) = {%nd,K\Sl\%nd,K) =0 



and 



(9a) 
(9b) 



where —m coincides with the average magnetization of the 
ground state. Consequently, the one-tangle has the form 



n{N, h) = l-4|(<I)gnd,^^|^^|$gnd,A0 



1- 



K[N, h] 



(10) 



In Fig. 131 Ti is shown as a function of /i/J, for different values 
of N . Similarly to the plots shown in Fig.I^P, it is worthwhile 
noticing that, for TV > 100, the curves are well approximated 
by the analytical formula 
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FIG. 3: One-tangle ri vs. h/J for several values of A'^. 
For N > 100, the plot is well fitted by the thermodynamic 
limit formula 
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obtained in the thermodynamic limit, and yielding the dominant behavior 

Ti{N oo,h^ h;) ^ (h, - h)^^^ . (12) 
B. Bipartite entanglement 

By taking into account the symmetries of the Hamiltonian, it is easy to show that the two spin reduced density 
matrix has the so called X-form and that the concurrence C{N, h, r) can be expressed in terms of the two-point 
correlation functions 

5^"^"'(iV,/i,r) = (ci>g„d,A'|5^^o5£,^J$g„,,A0 , (13) 

in which — x,y,z [l8l |. In particular, when two spins are at a distance r from each other, one gets C{N,h,r) = 
max{0, Ca{N, h, r)}, where 

CaiN,h,r) ^ A^/[g^^{N,h,r)]^ + [g^y{N,h,r)]^ - 2^/[l/A + g^'^{N,h,r)]^ - fh^ (14) 

denotes the antiparallel concurrence. Since the latter depends on the longitudinal spin-spin correlation function, 
g^^{N, h, r), its expression can be obtained analytically only in the range of h for which such a function is is known (see 
appendix and is otherwise computed numerically. 

For h ~ 0, using g^^{N,h = 0, r) = go{r), for each j — x,y,z, and g^y{N,0,r) = 0, one obtains the zero-field 
concurrence: 



C(iV, 0, r) =max{0,4|.go(r)| -2^1/4+ [<?o(0]'} • (15) 

An explicit evaluation shows that, for any N, Ca{N,h = 0, r) is always negative, except for the cases r — 1 and 
r = N — 1. The same property holds in the thermodynamic limit where go{r) takes the expression: [lo| 

9o{r) ~ —z — / dq^ dqp , (16) 

and the condition Ca{N, /i = 0, r) > is satisfied only for r = 1. Therefore, for h — 0, only the concurrence between 
spins lying onto nearest neighboring sites is different from zero, as it could be expected for an isotropic spin liquid 
state. 

Due to the lack of an explicit analytical formula for g'^^{N,h,r), the concurrence in the intermediate range < 
h < hs needs to be computed numerically. A result in the whole range can be easily found in the case < 20; 
in particular. Figs. |3]A. and Hp show the behavior of the concurrence vs. h, at fixed r, for A^ = 10 and A^ = 20, 
respectively. One clearly sees that the switching on of the magnetic field leads to an enlargement of the entanglement 
range, as C{N, h, r) becomes positive also for r > 1. 

For A^ > 50, it is still possible to determine C{N, h, r) exactly but only at values of h close to the saturation field. 
This is shown in plots of Figs. UP and|3)D, where the behavior of the concurrence near hg is shown for A^ = 50, 100 
and r = 1 — 5. In fact, when h ^ hg, the leading contribution to C{N,r,h) can be calculated by thermodynamic 
arguments, putting together the relationships (|A6|) and (|A11[) derived in the appendix For relevant values of 
r {r 7^ 0, A^), one obtains 



C,(N-,^,h,r)^4qo{hiqor)\~-J(l~^]T[qo,r]+(^] T^[qo,r]} , (17) 



where the functions 7 and are also defined in the appendix \X\ and go, the magnitude of the 'Fermi momentum' 
indicated in Fig. [l] can be expressed as go = ^(1 — 2m). As a result, one finds 



CaiN ^oo,h^ hg, r) ^ yJ{K - h)/J $(go, r) , (18) 

where <^{qQ,r) is a regular function of qa,r. In particular, at /i — ^ hs, one gets go = implying 7(0) = 7r/2 and 
T[qo — > 0, r] = 0, as it will emerge from the discussion in the appendix El Thus, 



a(iV^oo,/w/i.,r) = l^^^+o('^^^) , (19) 
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FIG. 4: C{N,h,r) vs. h at fixed r and for various values of A*'. Upon increasing h, C[N,h,r — 2) can be made > 0, thus 
yielding nonzero second neighbor concurrence. At a higher value of h, C{N, h,r > 3) becomes > 0, as well. 




FIG. 5: Behavior of C{N, r, h) for h hs with r = 1 (left panel) and r = 2 (right panel) and for N = 10, 20, 50, 100, and 
N ^ oo (thermodynamic limit). Upon increasing h, the entanglement range increases while at hs all the curves show the 
universal behavior of Eq. (|19|) 
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independently of r. The remarkable collapse of the curves onto each other for h hs, reported in Figl5]for r = 1, 2 
and for N = 10, 20, 50, 100, confirms the prediction of Eqs. (H?]), dH]). 

In general, one can say that by increasing the distance r, the concurrence C{N, h, r) is different from zero only in a 
smaller and smaller region of the magnetic field and that the range of the bipartite entanglement, rc, extends to the 
whole chain near the saturation point /i,. The behavior of rc vs h/ J is shown in Fig. [6] This is clearly reminiscent of 
what happens at factorization points 
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FIG. 6: Behavior of the range of the bipartite entanglement vs h/J; as h ^ hs, all the spinons are entangled and Vs — > N/2 



In accordance with this observation, we show below that near the saturation field, the bi-partite entanglement is 
dominant with respect to multi-partite ones, as it is the case for factorization points. 



C. Two-tangle and entanglement ratio 



In order to understand wether the entanglement stored in the ground state has a multi-partite rather than a bi- 
partite nature, we consider in this section the ratio between the two-tangle and the one-tangle defined above. Due to 
the monogamy relation for entanglement, this ratio is bound between zero (for purely multipartite correlations) and 
one (for bipartite entanglement) |17l |. 

We first focus on the two-tangle, that gives an overall measure of the bi-partite entanglement to which a given spin 
participate. It is defined by 

N/2 

T2{N,h)^2Y,C{N,h,r)^, (20) 

r=l 

and its behavior with respect to the magnetic field is reported in Fig. [7]A.. 
When h ^ hs, the limiting trend (black dashed line) is given by 

T2iN,h^h,)^{K-h)'', (21) 

with an exponent that we estimated, by fitting, to be a = 0.50 ± 0.02. This already implies that, near hs, T2 goes to 
zero as the one-tangle: n ^ {hs — ft.)"'^. 

To better show this fact, we finally consider the entanglement ratio R{N, h) = T2{N, h)/Ti{N, h), whose behavior is 
reported in Fig. [7p. Near the saturation field, the limiting trend R{N, h hs) ^ {h^ — h)'^ has an exponent (3 taking 
values smaller than 0.03, which is consistent with the other estimates given above. More importantly, the plot shows 
that at /i — >■ the one-tangle dominates, whereas in the neighborhood of hg the entanglement is essentially bi-partite. 



IV. CONCLUDING REMARKS 



Typical spin models considered so far to discuss the behavior of entanglement lie into two main categories: those 
with nearest neighbor interaction, such as the Ising model in a transverse field, and those with long range interaction 
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FIG. 7: A: Two-tangle, T2 as a function of the magnetic field. The black dashed line gives the limiting behavior reported 
in Ea. (|2ip . B: Entanglement ratio as a function of the external field. The plot suggests that near the saturation point, T2 
dominates over n. 



with spins residing on a fully connected graph, such as the LMG model or the uniaxial model [2l| . In the former 
case, bipartite entanglement decays quickly with the distance, and only nearest and next to nearest neighbor 
entanglement is found to survive. Notable exception are factorization points around which the concurrence is found 
to have an infinite range [l^. In collective spin models, on the other hand, the qubit-qubit entanglement is strictly 
zero in the thermodynamic limit, and its presence in the ground state is only due to finite size effects. In particular, 
the concurrence is found to scale as (TV — where is the number of spins. 

We have found that the inverse square (statistical) interaction among spins produces a quite different behavior 
of the entanglement as a function of the distance. In particular, we have shown that, in the absence of external 
magnetic field, h = 0, the system only displays nearest neighbor concurrence, while the entanglement has essentially 
a multi-partite nature; while, by increasing h, the entanglement becomes essentially bi-partite and the range of the 
concurrence increases. For a large enough h, the spin system eventually saturates, entering a fully polarized phase 
described by a completely separable ground state. We have discussed in particular the behavior of the entanglement 
near this saturation point, both for a finite size system and in the thermodynamic limit. In particular, we found 
that the divergence of the range of the bipartite entanglement is governed by the same exponent {\h — hs\^^^) of the 
ordinary isotropic Heisenberg model, consistently with the fact that the latter belongs to the same universality class 
as the Haldane-Shastry model. 

Since the inverse square interaction among spins in the HS model is essentially of statistical nature, the fact that the 
entanglement range tends to diverge near the saturation can be compared with the behavior of entanglement for other 
systems of indistinguishable free particles. As mentioned in the introductory section, the spin entanglement of a free 
electron gas is different from zero within a distance of the order of the Fermi wavelength. On the other hand, we have 
shown that in the presence of a large number of semions (i.e., near h ~ hs), the range of the bipartite entanglement 
extends to comprise the entire system. This behavior is ultimately due to the fractional statistics: by increasing the 
value of the external field, more and more spinons can be added to the system with smaller and smaller momenta, 
starting from the edge of the Brillouin zone and going towards its center. This implies that the only length scale of 
the system is the 'Fermi' wavelength q^^ , which, because of the semionic statistics, diverges when the saturation point 
is approached. In a sense, this behavior interpolates between the cases of entanglement among bosons and fermions. 
Free bosons can be always in a factorized state because there is no correlation between their relative distance and 
their polarization; for free fermions the range of bipartite entanglement between their spin is spatially limited to a 
finite Fermi length. For the semionic gas a Fermi length does exist, but it diverges at hg. Therefore, because of their 
quantum statistics an overlap of the spinon wave functions (whose extension is roughly given by h/qo) can take place 
at hs, ultimately making the saturation phenomenon a result of spinon condensation. 
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Appendix A: Magnetization and Spin correlation functions on \^gnd,K) 

This appendix contains the one- and two-spin correlation functions on j'l'gnd, ^"^')- These are the basic bricks to 
build the single- and two-spin density matrices needed to evaluate the entanglement. When possible, the calculations 
are carried out analytically on a finite chain and the thermodynamic limit is systematically worked out by sending 
N — > oo, while keeping K/N constant. When exact analytical results are lacking, namely in the calculation of the 
longitudinal two spin correlation functions, the corresponding quantities have been evaluated numerically. 



1. One-spin average values 

Since |$gnd,A') is an eigenstate of = J2a=i^a^ '^i^h total eigenvalue —K/2, one readily obtains 
(^gnd, ^I'S^I'i'gnd, ^''^) = {^gnd,K\Sy\^gnd,K) = 0. Morcovcr, bccause of the translational invariance on the lat- 
tice, one also has 

($g„d,if|5^|$g„d,/0 = -i + M J2 \'^k{1,z2,...,zm)\' = -^ . (Al) 

Z2,....Zm 

The thermodynamic limit is defined as N ^ oo, with K/2N ~ fh ~ constant. Therefore, for N — > oo, one can write 

{%^d,K\S'J%nd,K)=-fh . (A2) 



2. Transverse two-spin correlation functions 

The transverse two-spin correlation functions, introduced in Eq. (|13p . arc here denoted g^^ (r), for = x,y and 
r = l,...,iV — 1, where for notational simplicity the dependence upon N and h has been omitted. These functions 
can be expressed in terms of the correlation functions of the operators = 5^ ± iS*^ as 



(A3) 



with g\\{r) = g^^^^, K\Sq gi-id, K) . The auxiliary function g\\{r) can be exactly calculated by noticing that the 
act of flipping down one more spin in the state l^'gnd, K) is equivalent to the creation of a pair of 4,-spinons "on top 
of each other". Following the technique developed in Ref. 0, to compute the function g\\{r) for K = 0, one obtains 



M-l nia 
nia— m^— 

for a generic value of with 

_ [ma - mp + 1/2) T[M -ma,- \] r[m„ 1] r[m^-hl/2] T[M - mp] 



2N 



T[M- 



r[m„-h3/2] T[m0 + l] T[M -mf} + 1/2] 



(A4) 



(A5) 



To extract the asymptotic form of (r) in the thermodynamic limit, the Stirling's formula r[z] « ^/tt{z~1Y 2 e 
is used to approximate Xmo.ms ■ As a result, Eq. (jA4p becomes 



Mr 

167r 



go 



dqa dqp 



<?0 



\qa - qf}\e 



(A6) 



where qo is the 'Fermi momentum' introduced in Fig. [TJ An alternative formula for Eq. (|A6[) is obtained by using the 
integration variables (f>a and (/)/3, calculated from qa{p) — qosi'!i{4>a{i3))- In this way, one obtains 



with 



= 16^ 



Tr/2 



tt/2 



9\\{r) ~ {-lyqoliqar) , 



d(l)fi I sin(<?!)„) - sin((^/3)| e*P[''''^('^°)+"'"("^'')l 



(A7) 



(A8) 



Eq. (|A7I) is particularly useful for h ^ hs (that is, for go 0). 
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3. Longitudinal two-spin correlation functions 

The analytical computation of the longitudinal two-spin correlation functions, also defined in Eq. p3p for j = j' = z 
and r — 1, . . . , N — 1, is quite a formidable task not yet fully addressed [ij, [H,!!^. Even though in this paper g^^{r) 
is computed numerically for finite N,K, its leading contributions in goj ^-s well as its exact analytical expression for 
hs, are analytically determined. Indeed, by using the explicit form of \^gnd,K) given by Eq.([2]), one obtains 

g^^{r) = l-^{l-S,,,) + M(M-l) \'^K{l,e^^^^/'',Z3,...,ZM)\^ , (A9) 

where S'^ is set of the iV"' roots of the unity. Now, as pointed out in Ref [l5|, one can write 

J2 |<i>K(l,e2-/A^,Z3,...,^M)P = i^''-' n ^ :5^l<f>K(l,e2-'-/^,.3,...,ZM)P , (AlO) 

{z3,...,ZM}eS^ J=3 ^ 

for any integer L, in which 7 is the unit circle in the complex plane. By Eq. (jA10[) . one readily gets 

.9^^M = ^--(i-M+f— ) nao.r] . (All) 

4 TT \ / 

In this relationship, the function J^[qo,r] has the definition 

^ ^'^ ' {z3,....zj,/}GS2Af i<i=3 J=3 

for M > 2 (go > 0), and J"[go, = for M = 0, 1. It is straightforward to check that |J^[go,'^]| < 1 for any choice of 
go and r. Clearly, this gives a strong constraint on the behavior of g^^{r) for qq — >■ 0, that is, 

9''{r)^j-^{l-Sr.o)+0{ql) . (A13) 

4 TT 

Moreover, for small M, the function J^[qo,r] can be explicitly evaluated. For instance, for M = 2 {qo = 2tt/N), one 
gets 

.-w4-?(l-^.,o) + |^sin^(^) , (A14) 

while, for M ^1 {qo = n/N), the result is simply g^^i?-) ^ j - ^{1 - 5r,o)- 

In the absence of an external field, Q ^^ ( r) is exactly computed by simply using the fact that the ground state, 
I'J^gnd) if = 0), is a spin-0 spin singlet [l3,[3. Accordingly, one gets g^^ir) = g^^{r) = g^^{r), and g^^{r) = g^^^r) = 0. 
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